Controlling quasiparticle excitations in a trapped Bose-Einstein condensate 
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We describe an approach to quantum control of the quasiparticle excitations in a trapped Bose- 
Einstein condensate based on adiabatic and diabatic changes in the trap anisotropy. We describe 
our approach in the context of Landau-Zener transition at the avoided crossings in the quasi-particle 
excitation spectrum. We show that there can be population oscillation between different modes at 
the specific aspect ratios of the trapping potential at which the mode energies are almost degenerate. 
These effects may have implications in the expansion of an excited condensate as well as the dynamics 
of a moving condensate in an atomic wave guide with a varying width. 
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The trapped atomic Bose-Einstein condensate (BEC) 
provides a unique mesoscopic system in which the quan- 
tum mechanics of a gaseous system dominates. With the 
rich experimental techniques now available, we are at the 
point where "quantum engineering" of the condensate is 
an exciting frontier. In many examples of quantum en- 
gineering, strategically tailored states are engineered by 
controlling the population of the excited state of the sys- 
tem. For example, in cavity QED 0] this is achieved 
using the atom-held interaction to manipulate the pho- 
ton state, the atom state or a combination of both. On 
an atom chip 2] , time- varying or space- varying potentials 
are used to control atomic motional states. Turning back 
to BEC, it is well recognized that many-body quantum 
physics, and in particular the mean held theory, accu- 
rately describes most of the properties of the BEC. In 
this framework, many of the dynamical properties of the 
condensate can be understood in terms of the quasipar- 
ticle excitations of the system. The quasiparticle excita- 
tions also provide an effective tool for probing the role of 
interactions and for testing theoretical models. 

In this paper, we consider how changes in the shape of 
the harmonic trapping potential can be used to quantum 
control the condensate through the control of quasiparti- 
cle excitations. Such changes can be easily implemented, 
and may be encountered under different contexts such as 
the imperfections present in the guiding potential of an 
atomic wave guide for a moving condensate. Central to 
the our treatment is the control of the Landau-Zener (LZ) 
transitions between different types of quasiparticles. 
The LZ transition is observed in various systems, from 
atomic and molecular physics 0, to solid state physics 
through to mesoscopic systems 0. However, in general, 
the theory only applies well to linear systems. By con- 
trast, the successful application of LZ theory to a BEC 6] 
requires more care given the intrinsic nonlinearity of the 
mean field that describes the condensate. 

We begin with a brief description of the hydrodynamic 
modes in an anisotropic trapping potential based on the 
observed structure of the quasiparticles. We then exam- 
ine the LZ transition qualitatively and quantitatively. We 
additionally find that population oscillations between dif- 
ferent types of excitations around the avoided level cross- 



ing is possible. Methods to control the normal modes are 
discussed. 

To describe the BEC ground state, we consider a 
pancake-shape trap geometry with a strong confinement 
along the z-axis Q- A two-dimensional (2D) time- 
independent Gross-Pitaevskii equation (GPE) Q based 
on mean field theory is used: 

(H + 9^0-^)^0 = 0, (1) 

where H = -(h 2 /2M)\7 2 + V tT {r). Here, ^ (r) = (^(r)) 
is defined following the standard mean field treatment, 
with i}}(y) being the field operator satisfying the bosonic 
commutation relation, g = Anah 2 /M is the interparti- 
cle coupling constant with the s-wave scattering length 
a. fi is the chemical potential, and the wavefunction 
■00 is normalized to the total number of particles N. 
For our calculation, 2D anisotropic harmonic potential 
Vtr = M(uj 2 x 2 + uj 2 y 2 )/2 is adopted and an effective 
.92D = g(mu} z /2irh) 1 / 2 is used Q. Coupled Bogoliubov- 
de Gennes (BdG) equations have been solved numer- 
ically to calculate low-lying energy spectrum at various 
aspect ratios, A = Lj y /uj x : 

(H + 2gipQ - m) u j - 9^o v 3 = tWjUj 
{H + 2g%yb - /z) Vj - g% 2 u 3 = -HujjVj, 

where Uj, Vj and u>j represent the wavefunctions and 
eigenenergies for the jth quasiparticle excitation. Time- 
dependent GPE [8j has been used for the examination of 
dynamics with varying parameters. All values in this pa- 
per are based on the dimensions of trap potential along 
the x axis: yjTi/Muj x for the length, and lu^ 1 for time. 

The quasiparticle energy spectrum from solving BdG 
equations for different values of the aspect ratio, 1 < A < 
5, is shown in Fig. |T] As A increases, the BEC under- 
goes an effective dimensional change from 2D to almost 
ID. On the left end of the plot, the aspect ratio of the 
trap is 1, for which the full spectrum had been found 
[lfj. In this case, every energy level has a degeneracy of 
2 and is characterized by the angular quantum numbers 
±m and a radial quantum number n. If the aspect ratio 
deviates from unity, this degeneracy is lifted: the degen- 
erate states (n, ±m) become two non-degenerate states 
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FIG. 1: Energy spectrum of hydro-dynamic modes of BEC 
with respect to the aspect ratio X = ui y /u x . In the circle is 
the avoided crossing between (10,0) and (2, 2) modes. 

(n x = \m\ — l,n y = 2n + 1) and (n x = \m\,n y = 2n). 
The new set of quantum numbers (n x ,n y ) corresponds 
to the number of zeroes along each axis in the quasipar- 
ticle wave functions u.j and Vj . Examination of the time- 
dependent probability density \tpo + Uje~ W3t - v*e luJit \ 2 , 
where ipo indicates the ground state wave function, shows 
a standing wave with n x and n y nodes along each axis. 
In Fig. ^ the quantum numbers (n x ,n y ) are indicated 
for each line, and it is noted that different slopes corre- 
spond to different n y , as we are using the u) x trap dimen- 
sion. A number of avoided crossings between spectral 
lines are clearly observed in Fig. ^ Since the "pertur- 
bation" added to the Hamiltonian due to the aspect ra- 
tio change, M(x 2 u) x 5u x +y 2 uj y 5uj y ), involves even parity, 
these avoided crossings are found to occur only when two 
non-crossing states have the same parity llj. It should 
be noted that avoided crossings are found regardless of 
the dimensionality of the system [l^. For high A, in the 
low energy regime, n y = modes are decoupled from 
the others to form a group of ID excitations, for which 
the wavelength of the excitation along the x axis is much 
larger than the size of BEC along the y axis. 

We focus our attention on a specific avoided crossing 
between states (n x = 2, n y = 2) and (n x = 10, n y = 0) 
which is circled in Fig. 2] Figure [21 shows the change 
of the function n'j = ip*Uj — t/jVj Q for an energy eigen- 
states along the lower spectral line of the avoided cross- 
ing. Figure|21(a) and (c) correspond to (2,2) and (10,0) 
modes with A = 3.1 and 3.5 respectively. We shall re- 
fer to the asymptotic states far from the center of the 
avoided crossing such as (a) or (c) "shape eigenstates" as 
they have definite number of zeros along each axis. The 
energy eigenstate at the avoided crossing, Fig. |21 (b) , is a 
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FIG. 2: Change in the density fluctuation n' for the energy 
eigenstates along the lower spectral line associated with the 
avoided crossing in the circle drawn into Fig. 

linear combination of the two shape eigenstates, (a) and 
(c), and is itself not a shape eigenstate. 

If the system undergoes an adiabatic change in its 
aspect ratio around the avoided crossing starting from 
one of the asymptotic states, it cannot go through the 
energy gap between the two avoided spectral lines 0, 
resulting in the exchange of the mode characteristics; 
for instance, starting from a (2, 2) mode, it ends up 
in (10,0) mode. For a diabatic transition, the state 
may tunnel through the energy gap between the two 
avoided spectral lines to remain in the same shape eigen- 
state. According to the LZ theory the probability 
for the state to remain in the same shape eigenstate is 
P = e- 2 ^, where 7 = TT l {&/2) 2 [{d/dt){e x - e 2 )] _1 . 
Here, ei and e 2 are the eigenenergies corresponding to 
two avoided states and A is the minimum value for 
| ei — £2 1 . For the above avoided crossing A w 0.07 and 
(d/dt)(ei - e 2 ) = \{djd\){e\ - e 2 ) ~ 1.4A, which gives 
P = e - 5 - 5x10 A in the u> x trap dimension. Here, A 
denotes dX/dt and A = 0.001 corresponds to the adia- 
batic case with P = 0.004, while A = 0.1 to the dia- 
batic case with P — 0.95. In order to observe this adia- 
batic transition experimentally, the change in A from 3.1 
to 3.5 should happen throughout the time longer than 
At = AA/A = 400 w 64 x 27TW" 1 . 

We have numerically simulated the temporal evolu- 
tion of a chosen excitation with A varying in time by 
using time-dependent GPE. As an initial state, we chose 
the normal mode (2,2) for A = 3.1. To visualize more 
clearly the excitation possessed by the condensate as it 
evolves, we subtracted the ground state from the total 
wavefunction to give time-dependent density fluctuation 
8p(t) = ip*(t)ip(t) — tpo^o- Assuming that there are 
\qj\ 2 quasiparticles for each excitation mode j, as far as 
\qj\ 2 <C N, the time evolution of the mean field is almost 
exactly given by [3 

m ^o + £ (we-**** ~ q*v*^ 1 ) , (3) 

3 

which, with real Uj, Vj, ipo an d qj, yields 

*p(t)»5^2^nJco8(« J t), (4) 

3 

with higher order terms of qj ignored. Noting that each 
term in 5p(t) is a product of n'j with a sinusoidal time 
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FIG. 3: Distribution function in the frequency space for each 
value of A while A increases from 3.1 to 3.5 (a) adiabatically 
(A = 0.001), and (c) diabatically (A = 0.1). (b) is for the 
intermediate case (A = 0.01). 



dependence cos(cjj-t), we have compared Sp(t) with the 
exact energy eigenstates given in Fig. [21 We have con- 
firmed that the shape of Sp(t) changes from Fig. 0a) to 
(c) under the adiabatic change in A, while the state re- 
mains in the same shape eigenstate for the diabatic case. 

In order to estimate the amount of transition quantita- 
tively, we have used spectral analysis. Equation shows 
that the quasiparticle population for each energy mode 
can be obtained by performing a time-domain Fourier 
transformation of Sp(t). As A changes in time the quasi- 
particle populations change as a function of A. Figure 
[3] shows the distribution function in the frequency space 
for each A. For the adiabatic case (a), the peaks nicely 
follow a single spectral line of two avoided energy levels, 
while for the diabatic case (c), the population peak tun- 
nels through from one spectral line into the other to stay 
in the same shape eigenstate. In addition, we obtained 
the result for the intermediate case (b) at A = 0.01, which 
shows the splitting of the initial population on one en- 
ergy eigenstate into two after going through the avoided 
crossing region. It is notable that, although these sim- 
ulations were carried out using full GPE, the linear LZ 
theory is still valid with the quasiparticle populations 
small enough compared to the ground state population. 

We believe that the rate dependence of the process may 
have an implication in the process of releasing conden- 
sates with specific excitations from a trapping potential, 
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FIG. 4: (a) The population growth of the (10,0) mode while 
the ground state of BEC is driven by Vd T (t) ~ sin(wt). The 
other mode (2, 2) is clearly suppressed, and (b) the follow-up 
population oscillation between (10, 0) and (2, 2) modes near 
the avoided level crossing. 



which is usually done in the experiments. In addition, 
the roughness of the potential in the atomic waveguide 
could introduce the change of aspect ratio naturally, in 
which case, the reaction of the phonon modes may have 
a strong dependence on the velocity of a moving conden- 
sate. 

We now consider an additional scenario for quantum 
control. With a correct driving potential it is possible 
to excite energy eigenstates or certain shape eigenstates. 
However, if the generated shape eigenstate coincides with 
an avoided crossing for some specific aspect ratio A, the 
excitation is not an energy eigenstate but a linear com- 
bination of two almost-degenerate energy eigenstates. If 
one denotes the two energy eigenstates as Xi an d the two 
orthogonal shape eigenstates as <pi, these two different ba- 
sis sets are related by a unitary transformation Uij . The 
temporal evolution of an initial shape eigenstate, ip(t = 
0) = fa = Unxi + U 12X 2, is then ^{t) = UnXie^* + 
Ui2X2e~ luJ2t , such that the population of mode 4>i is 
|(^i|V(i))| 2 = |C/n| 4 + |;7 12 | 4 + 2|C/ 11 | 2 |C/ 12 | 2 cosK-c l ; 2 )i, 
while \{(fa\ip(t))\ 2 = 1 - \{(j>i\ip(t))\ 2 . Clearly the popula- 
tions oscillate between two modes even though it starts 
from only one of the two states. In our numerical simu- 
lations, we could efficiently generate a shape eigenstate 



with a specific quantum number (n x 



by applying a 



driving external potential Vd r (t) to the trap potential Vt r 
of the form: 



V Al (t) 



. . . f2irx\ f2ny\ 
smM)cos^— jcosj^— j, 



(5) 



where u), X x and X y are the frequency and the estimated 
wavelengths for the aimed normal mode respectively. A 
(10,0) mode at A = 3.3 was generated using Vd T (t) with 



6.2, X x w 3 and A^ 



oo. When Vdr was turned 



off, the oscillation between the generated state and the 
(2, 2) mode was clearly observed in the evolving density 
fluctuation p(t). In order to show this "shape oscilla- 
tion" quantitatively, we adopted quasiparticle projection 
method [13j, in which the amplitude of the jt\i mode 
from Eq.~ © is qj = J dr(u*^(t)e 1 ^ + w|V*(*)e~ iM *), 
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FIG. 5: Schematic diagram of an avoided crossing. <f>i indicate 
asymptotic shape eigenstates and \i exact energy eigenstates. 

where iij = Uj — ajipo/N , v* = v* — a* V^o/Ao and 
dj = J dv^Uj = J dripoVj. Fig. ^a) shows the pop- 
ulation growth of the (10, 0) mode in the BEC driven by 
Vdr(i), while the population of the (2,2) mode is sup- 
pressed. Fig. 0Jb) represents the follow-up population 
oscillation between two shape eigenstates after turning 
off the driving potential. The period of the shape oscil- 
lation matches with the expected 2ir/(uii — 0J2) ~ 90. A 
related result was reported in Ref. [13 which was dis- 
cussed in a different context of resonant nonlinear mode 
coupling. It is noted that a similar oscillation can also be 
found in particle physics, in which neutrinos, generated 
in a muon or electron eigenstate, are slightly rotated from 
the exact energy eigenstate and undergo what is known 
as "neutrino oscillation" [T^| . 

Finally, we would like to mention further possible 
quantum control of BEC involving phonon modes around 
avoided crossing. We have observed that the transition 
between two different shape eigenstates happens under 
adiabatic process. The same kind of transition, how- 
ever, may also be made with a different process (See 



Fig. [SJ|: Starting with a shape eigenstate 4>i at A = Ai, 
make a diabatic or sudden change of A to A ac which does 
not change the mode, wait until the shape oscillation 
achieves </>2 > and then make a diabatic change of A again 
from A ac to A2 ending up with the other shape eigen- 
state cf>2 at A2 (G — » A — > D). If we do the second 
diabatic process backward from A = A ac to Ai, we are 
able to get mode transition from cf>i to 02 at the same 
A = Ai only by changing the aspect ratio slightly back 
and forth (G — > A — > E). Another possible process is 
to change A from Ai to A2 and back to Ai with a fixed 
|A| and with the shape eigenstate (f>\ initially. The fi- 
nal state will end up with the same shape eigenstate for 
both diabatic (G^A^F^A^G) and adiabatic 
(G^B^D^B^G) case, except for the interme- 
diate case, in which the final state will have considerable 
population on both shape eigenstates. This process may 
be compared to the collision of the second kind in an 
atomic collision [llJ. These could be useful if the BEC 
phonon modes could ever be used in the future as qubits 
for quantum computation. 

In summary, we have studied ways to quantum control 
quasiparticles in BEC with the aspect ratio of the har- 
monic trapping potential as a parameter. Using numeri- 
cal simulations, we confirmed that LZ transition does oc- 
cur in the linearized regime. We have also observed that a 
shape oscillation between different types of collective ex- 
citations are generated at a specific trap anisotropy. This 
"beating" observed in the condensate density may also 
be viewed as a signature for the existence of an avoided 
crossing in the Bogoliubov spectrum. These phenomena 
should be readily observable experimentally. Future work 
would involve extension of this work to a rotating BEC 
with vortex lattices |ln| . This work is supported by the 
NSF, DOE, ONR and ARO. 
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